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Floquet topological insulators are noninteracting quantum systems that, when driven by a time- 
periodic field, are described by effective Hamiltonians whose bands carry nontrivial topological 
invariants. A longstanding question concerns the possibility of selectively populating one of these 
effective bands, thereby maximizing the system’s resemblance to a static topological insulator. We 
study such Floquet systems coupled to a zero-temperature thermal reservoir that provides dissipa¬ 
tion. We find that the resulting electronic steady states are generically characterized by a finite 
density of excitations above the effective ground state, even when the driving has a small amplitude 
and/or large frequency. We discuss the role of reservoir engineering in mitigating this problem. 


I. INTRODUCTION 

In recent years, the possibility of engineering topo¬ 
logical states of matter in otherwise trivial materials 
has motivated significant interest in electronic systems 
driven periodically in time m- The prescription for 
a particular target topological state is obtained using 
Floquet theory mm, which describes a periodically- 
driven quantum system with Hamiltonian H{t) = 
t) in terms of a time-independent Hamiltonian Hp. 
In particular, Floquet’s theorem states that there ex¬ 
ists a complete basis of solutions of the time-dependent 
Schrodinger equation H{t)\ipj{t)) = idtlifjit)) of the 
form 

= \u,{t)), (1.1) 

where the state \uj{t)) = \uj{t + t)) is periodic in time 
with the same period as H(t). The quasi-energies Cj 
can be thought of as the energy spectrum of an effective 
“Floquet Hamiltonian” Hp defined by evaluating the 
evolution operator U (^ 2 , H) at an integer multiple of the 
period r, i.e. U{nT,0) =: exp(—i i/p ^t). Floquet’s the¬ 
orem guarantees the existence of a time-periodic unitary 
operator P{t) that maps the time-dependent Hamilto¬ 
nian onto Hp [MU]: 

Hp = PHt)H{t)P{t) -iP\t)dtP{t). (1.2) 

If the single-particle Hamiltonian Hp and its associ¬ 
ated eigenstates \uj{nT)) define a model with nontrivial 
Berry curvature, then it is possible for the quasienergy 
“bands” to possess nonzero Chern numbers or Z 2 in¬ 
variants [HE], defined by analogy with static Hamilto¬ 
nians m- 

Despite the formal similarity to undriven systems, 
there are several fundamental differences between a Flo¬ 
quet effective Hamiltonian Hp and a static Hamilto¬ 
nian. First, there is an ambiguity in the definition of 
Hp above: the quasienergies Cj are only defined modulo 
multiples of the driving frequency to = 2'k/t. In other 



FIG. 1: Ambiguity in the ordering of quasienergy bands. For 
two quasienergy bands falling within a given strip of size w, 
as in a), a gauge transformation of the form | [1.3[ ) can be 
used to shift the window and invert the ordering, as in b). 

words, the definition of Hp is invariant under the gauge 
transformation 

Cj —>■ Cj + rij (jj and \uj{t)) —>■ \uj{t)) , (1.3) 

for any set of integers {rij}. This makes an a priori def¬ 
inition of a “lowest” quasienergy impossible, since one 
can always fold and/or reorder (c.f. Fig.[^ the spectrum 
by means of such a gauge transformation m- 

This obstruction to defining a unique Floquet ground 
state highlights a second important contrast with static 
systems, namely the fact that there is no universal prin¬ 
ciple determining the occupations of quasienergy states 
at zero temperature. Indeed, this signature of the inher¬ 
ently out-of-equilibrium nature of periodically-driven 
systems poses a challenge to theoretical proposals of 
Floquet topological states — even if a Floquet system 
has topological quasienergy bands, there is no guarantee 
that the long-time state of the system is a pure Floquet 
state with the desired features [smiiiiH]- For open 
systems in contact with a thermal reservoir, as is the 
case in solid-state realizations, the detailed properties 
of the reservoir and its coupling to the system play cru¬ 
cial roles in determining the nonequilibrium steady state 
of the system at long times [TU . 

In this work, we investigate the occupations of topo¬ 
logical Floquet bands in a noninteracting fermionic sys- 
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tern coupled to a thermal reservoir at zero temperature. 
We ask whether, and under what conditions, it is pos¬ 
sible for the nonequilibrium steady state of the system 
to feature a single fully populated Floquet band when 
the original undriven system is half-filled. Using a Flo¬ 
quet master equation approach, we examine as a func¬ 
tion of the driving amplitude and frequency the rate 
to escape a given Floquet state. We find that generic 
systems undergo heating that spoils the complete oc¬ 
cupation of a single band, even in the favorable limit 
of weak, off-resonant driving. In the best-case scenario, 
the excitation density falls off as a power law in 1 /w, 
with the non-universal decay exponent set by the low- 
energy behavior of the bath density of states. We close 
by suggesting bath engineering schemes that could fur¬ 
ther suppress excitations. While our motivation derives 
from the pursuit of topological phases in driven systems, 
we point out that our results hold equally well in “triv¬ 
ial” Floquet band insulators. 

The structure of the paper is as follows. In Sec. |n) 
we provide an overview of the Floquet master equation 
approach used in this work. Rather than focusing on ex¬ 
position, this section provides a fresh perspective on the 
subject by emphasizing the requirement of invariance of 
all physical quantities under gauge transformations of 
the form (1.3). This physical principle resolves many of 
the conceptual ambiguities regarding the definitions of 
Floquet states and the associated quasienergies. Within 
this formalism, we outline conditions under which open 
Floquet systems can reach steady states that resemble 
their equilibrium counterparts, i.e. where the occupa¬ 
tions of the Floquet states at finite temperature are dis¬ 
tributed in a Boltzmann-like fashion. In the remaining 
sections, which can be read independently of Sec. [llj 
we specialize to the study of noninteracting, two-band 
Floquet systems coupled to zero-temperature bosonic 
reservoirs, and present the main results of this work, 
which were outlined above. 


II. FLOQUET MASTER EQUATIONS 

A. Definitions and comments on gauge invariance 

In this paper, we will study Floquet systems that are 
described by Hamiltonians of the form 

H{t)=Hs{t) + HsB+HB, ( 2 . 1 ) 

where Hs{t) — Hs{t + t) describes the periodically 
driven system, Hb describes the reservoir (or “bath”), 
and Hsb describes the coupling between them. Follow¬ 
ing Refs. Uni and 1231 we assume a system-bath coupling 
of the factorized form 


where 7 is a real coupling constant [assumed to be 
smaller than any energy scale in and where S 

and B are Hermitian operators acting solely on the de¬ 
grees of freedom of the system and the bath, respec¬ 
tively [25] . System-bath couplings of this form are ubiq¬ 
uitous in the study of open quantum systems, e.g. in 
spin-boson-type models |2Z]. The reservoir described 
by Hb, whose eigenstates and energies Hb |j^) = E,,\iy) 
are known, is assumed to be in thermal equilibrium (ei¬ 
ther at zero temperature or at a finite temperature 1 //?, 
although in Sec. |III| we will specialize to the case of zero 
temperature). 

The unitary time evolution of the full closed system 
is completely characterized by the density matrix g{t), 
whose equation of motion is given by 

idtgit) = [H{t),g{t)]. (2.3) 

Our interest, however, is primarily in the influence of 
the reservoir on the system described by Hs{t). The 
important quantity to consider is then the reduced den¬ 
sity matrix 


p{t) := Tre [g{t)] 


(2.4) 


where Tre) ■ ] represents the operation of taking the 
trace over the bath degrees of freedom. The evolu¬ 
tion of the open system described by p{t) is obtained 
by “tracing out” the bath in Eq. ( |2.3[ ), which results in 
a master equation for the reduced density matrix. Ow¬ 
ing to the explicit time-dependence of Hs, this master 
equation has time-periodic coefficients, so that the Flo¬ 
quet states defined by Eq. (l.I) form a natural basis in 
which to resolve p{t) naiiM]. The derivation of the 
master equation for the reduced density matrix in the 
Floquet basis proceeds via the Born-Markov approxima¬ 
tion, which amounts to an assumption of weak coupling 
between the system and bath, as well as the assumption 
that the correlation time of the bath is sufficiently short 
that memory effects can be neglected. 

A particularly clear derivation of this Floquet master 
equation, as well as a careful discussion of the hierarchy 
of time-scales necessary in order for the Born-Markov 
approximation to hold, is presented by Hone et al. in 
Ref. 1231 We refer the reader to that work (and references 
therein) for details, and begin our discussion with the 
master equation written in their notation: 

{dt + iezj) Pij{t) = ~2 ^ ' {pij{t)Eik;lk{t) + Pil{l-)Rjk-,lk{^) 

k,l 


— Pklit) [Rlj-ki{t) + Rki-lj{^)\ } ) 


(2.5a) 


where 


Py (t) := {Ui{t)\p{t) \uj{t)) , 


Hsb=7SB, 


( 2 . 2 ) 


iij ■— Si Sj, 


(2.5b) 

(2.5c) 
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and the coefficients Rij-ki{t) '■= Sk 
where 

Rij-,ki{K) := 2Tr^^y^Sij{m + K) Sli{m) g{eik - muj) 


Sij{m) ■■=- [ 

Jo 


{u,it)\S\u^it)). (2.5d) 


The influence of the bath is contained in the function 
g{E), essentially a weighted density of states, which is 
defined by 


9{E) = 6iE + E, - E^), 


Z 


fljl' 


^ = E< 


,-/3E. 


(2.5e) 

(2.5f) 


are the time-averaged rates. Equation (2.9a) can be 


further simplified by considering the structure of the 
averaged density matrix pij. Indeed, so long as = 
Ci — ej 0 for i ^ j, the system-bath coupling 7 can 
be chosen to be much smaller than the smallest eij. In 
this case, the off-diagonal elements of the steady-state 
reduced density matrix vanish to order 7^ (c.f. [23]), 
and the occupations Pi := pu of the Floquet states are 
solutions to the rate equation 

0 = [Rj^i Pj - Ri^j Pi), (2-10) 

where the transition rates 

Ri^j = Rij;ij = Bij-ij 

= g{eji-muj). ( 2 . 11 ) 


where j3 is the inverse temperature. 

Equation (2.5a) is advantageous in that it is com¬ 


pletely invariant under the gauge transformations (1.3) 


that reshuffle the quasienergy spectrum. Indeed, defin¬ 
ing 




S := ej + rijUJ, 

one finds that the density matrix transforms as 


(2.6a) 

(2.6b) 


Consequently, the left- and right-hand sides of Eq. ( |2.5a ) 
transform with an oscillating phase factor that cancels, 
and the equation is therefore invariant under the gauge 
transformation. 

Despite the appeal of a gauge-invariant equation of 
motion for the reduced density matrix, it is also difficult 
to make analytical progress while the rates Rij;ki{t) are 
time-dependent. If one is only interested, as we are, 
in the long-time limit, where the density operator does 
not vary substiantially over a single period, then a valid 
means of bypassing this diffic ulty is to consider the time- 
average of both sides of Eq. (2.5a) [33], 


Pij — „ 'y ' [pijRik;lk + PilRjk-,lk 


k,i (2.9a) 

~Pkl {Rlj-,ki + J 


where pij = Pij(t —)■ 00) is the time-average of Pij(t), 
and where 


Rij-kl — Rij-,kl(^E — 0) 


(2.9b) 


Our treatment of nonthermal steady states in Sec. m 
will be based on a careful analysis of the rate equation 
(2.10) in the special case of a two-state system. 


The time-averaging p rocess outlined above spoils the 
gauge-invariance of Eq. (2.5a); after time-averaging, one 


therefore implicitly commits to a choice of ordering for 
the quasienergies. Indeed, observe that 


Rij;kli_E^ — Rij-jkl^E -f riik 


( 2 . 12 ) 


p,^.(t) = pC(t) = (umPit) \u',it)) , n,k - 

(2.7) 

while the time-dependent rates transform as 

= ( 2 . 8 ) 


In particular, Rij-ki{K = 0) ^ B'ij-,ki{E = 0) if 
riji ^ 0. This indicates that the primed ver¬ 


sion of Eq. (2.9a) need not hold given the unprimed 
version. However, observe that the diagonal rates 
Rij-^ij = Ri^j, along with the diagonal density matrix 
elements Pii{t), are invariant under the gauge transfor¬ 
mation (2.6) [c.f. Eq. (2.7)]. Therefore, the rate equa¬ 


tion (2.10) that governs the steady-state populations is 


gauge invariant. 

One important consequence of the gauge invariance 
of Eq. (2.10) is that, in addition to exact degenera- 
one must also treat carefully quasi- 


cies where 


= 0 , 


degeneracies, where e^- = mu for some m G Z. Such sit¬ 
uations were of primary concern in Ref. 1231 which points 
out that such degeneracies can have profound effects 
on the steady-state density matrix. In the language of 
gauge-invariance adopted here, such quasi-degeneracies 
are simply indicative of the fact that there exists a gauge 
in which quasienergies and are degenerate. 


Consequently, in such situations, Eq. ( [2.10] ) is no 
longer valid in the basis of Floquet states. However, the 
density matrix nevertheless has a block-diagonal struc¬ 
ture, Pij being nonzero only when tij = 0 in some gauge. 
(Let us, for simplicity, adopt this gauge here and for the 
remainder of the present discussion.) To close this sec¬ 


tion, we show that a rate equation of the form (2.10) is 
recovered upon transforming to the basis in which pij 
is diagonal. The unitary transformation M that diago¬ 
nalizes the density matrix acts on the Floquet states as 
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B. Conditions for thermal Floquet steady states 


^ \n{t)) . 

i 

p., = E ni. 

s« = E 




Rirki= V M,tM^rR-. 

^ It 33 kk II 13] 

i,j ,k,l 


kl 


(2.13a) 

(2.13b) 

(2.13c) 

(2.13d) 


We now plug these transformed quantities into 
Eq. (2.9a), keeping in mind that its left-hand side can 


be set to zero given that pij is zero whenever ^ mw. 
We obtain, for each i and j, 


0 = E («a j) E [m ihm + h 

ij k,l 


hi {hj-M + ^h-ij) 


(2.14) 


One can show that Eq. (2.14) holds if and only if 


® “ X! [h hk-ik + hi ^]k-ik hi (%;fci + 

k,l 

. 

for all i and j. Next, consider the above equation when 
* = J) while at the same time keeping in mind that p is 
diagonal. We obtain 


^ — '^^hi {kik-~ik + ^Ik-^k 


or equivalently 


"X! hk ij^ki-ki + ^ki-ki 


^^J 2 {h^h^k-Pkh^^l - 

k 


(2.16) 

(2.17) 


where pi is the probability of being in the state labeled 
by i, and 

Ri^k = ‘2T^hJ2\hki'rn)\^ 9hki-^)- ( 2 - 18 ) 

m 

This discussion has therefore demonstrated that a rate 


equation of the form (2.10) always determines the 


steady-state reduced density matrix pij in some basis, 
namely the basis in which pij is diagonal. Whether or 
not this basis is the basis of Floquet states, or some 
other time-periodic basis, depends on the degeneracy 
structure of the quasienergy spectrum. Nevertheless, 
the subsequent analyses presented here hold equally well 
in this choice of basis. 


In this section, we will determine, within the master 
equation formalism discussed above, conditions under 
which the driven system relaxes to a thermal distribu¬ 
tion with respect to the quasienergies. In particular, 
we will show that such a situation occurs provided that 
pl(t)i7sB-P(I), where P{t) is defined via Eq. (1.2) with 
H{t) = Hs{t), is either time-independent or depends 
on time in a particularly simple way. While statements 
to this effect have been made in various works [25H3T], 
we provide here a simple and complementary derivation 
of the statement from the master equation formalism 
outlined in Ref. [23) . and provide connections to the 
principle of gauge invariance discussed previously. Be¬ 
fore beginning with the derivation, we observe that, at 
finite temperature, the function g{E)^ which contains 
the influence of the bath on the driven system, satisfies 


g{E)=gi-E)e-f^^, 


(2.19) 


which induces the following generalized detailed balance 
relation for the rates appearing in Eq. (2.10) [23]: 


R 


R, 








E urn p-0muj ’ 
m 


where 


R7hi = 27r7^ \Szj{m)\ g{ej, - mw). 


(2.20a) 


(2.20b) 


Note that substituting Eq. (2.20) into the rate equation 


(2.10) yields a thermal distribution for the occupations 
of the Floquet states, i.e. pj oc if m = 0. 

We can now show that such a thermal distribution 
emerges if Pl'(t)i7sB-P(0 is time-independent, simply 
by showing that it implies that only the m = 0 terms 
contribute to Eq. ( |2.20 ). To do this, we make use of an 
explicit representation of the operator P(f), which can 
be derived as follows. Noting first that the evolution 
operator U (^ 2 , ^i) can be written in terms of the Floquet 
states \uj{t)) as 


U{t2,ti) = ^ |V'i(t2)) 

3 


( 2 . 21 ) 


we compute 

idtU{t,Q) ^ Hs{t)U{t,Q) 


= E' 


ej \uj{t)) (uj(O)l 


+ idt\uj{t)) {uj{0)\ . 


( 2 . 22 ) 
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Acting from the right with W{t,0) on both lines above 
and using the fact that {ui{t)\uj{t)) = Sij at any time 
t, we deduce that 

Hs{t) = [^3 ■ 

3 

(2.23) 

At this point, making the ansatz 

P(t):=^|u,(t))(n,(0)|, (2.24) 

3 

we find that indeed 

Ht? = P\t)Hs{t)P{t) - iP\t)dtP{t), (2.25a) 

where 

(2-25b) 

3 


as desired. 


Using the explicit form of P{t) provided in Eq. (2.24), 
the desired result follows directly from the definition of 
given in Eq. (2.20b). In particular, observe that 

S^J(.m) = - {u3t)\S\u,{t)) 

T Jo 


'M0)\PHt)SP{t)\u,{Q)). 

(2.26) 




dt e 


— imcot 


If P^t)SP{t) [and therefore P^(t)i?sB-P(t)] is indepen¬ 
dent of time, then Sij{m) = = 0 for all m 7 ^ 0. We 

have therefore shown that, if the operator P^t)HsBP{t) 
is time-independent, then the steady-state occupations 
of the Eloquet states are given by pj oc e~^'^3 ^ as they 
would be for a system with Hamiltonian Hp at equi¬ 
librium with a finite-temperature reservoir. The zero- 
temperature limit of the associated grand-canonical dis¬ 
tribution determines an unambiguous ordering of the 
quasienergies. 

If P^ {t)SP{t) is not time-independent, it is still possi¬ 
ble to reach an effective thermal distribution if detailed 
balance, namely the relation 


Pi->j _ % 
Pj^i Pi 


(2.27) 



FIG. 2: (Color online) Closed loop in a graphical represen¬ 
tation of the state space defined by the Floquet states and 
the rates Ri-^j connecting them. Blue circles represent Flo¬ 
quet states, and a directed arrow pointing from state i to 
state j represents the rate Ri^j. If the condition ( |2.31[ ) is 
satished around all closed loops in the graph, then the extra 
factors of mijui on the li nks can be removed by a gauge trans¬ 
formation of the form ( 2 . 61 , yielding a thermal quasienergy 
distribution. 


Since S' is a Hermitian operator, we have rriij = —m-a , 
and Sij{mij) = S*3mji). It follows that Eq. (2.20) 
becomes 


Ri 


R 


^^3 _ ^ 


(2.29) 


In this case, one recovers a thermal quasienergy distri¬ 
bution only if it is possible to consistently absorb the 
extra factor of rriijUj into a redefinition of the quasiener¬ 
gies ej. In principle, this can be achieved by means 
of gauge transformations of the form ( 2 . 6 ), which shift 


the quasienergies by integer multiples of the driving fre¬ 
quency Lo. In practice, however, it might arise that 
such a gauge transformation cannot be performed con¬ 
sistently over all Floquet states. Fortunately, there is 
a condition, demonstrated below, on the integers rriij 
that guarantees that such a transformation can be car¬ 
ried out. 

The set of nonvanishing transition rates Ri^j defines 
the directed edges of a graph whose vertices are the Flo¬ 
quet states (see Fig. [^. In order to consistently gauge 
away the extra factor of w, it is sufficient to require 
that for any closed loop —>■ 12 —t zs —>■••• —t ik-i —t 
ik —t *1 in this graph, one has 


holds for all i and j. For example, let us expand the 
time-periodic operator P^t)SP{t) and 

consider the case where (ni( 0 )| |'Uj( 0 )), for fixed i 

and j, is nonzero for a single mode n = rriij. In this 
case, 


R^^j = 2717^ \S^j{mij)\'^ g{€ji - rrUjUj). (2.28) 


Rii 


>Z2 Ri2^i3 


Ri2^il Ri3^i2 


Rik-l^ik Rik—iil 
Rik^ik-l Rii^ik 


= 1 . 


(2.30) 


[Observe that this condit ion holds automati cally if de¬ 
tailed balance, Eq. (2.27), holds.] Using Eq. (2.29), one 
finds that this condition holds if and only if 

miii2 + H-1- mik-iik + (2-31) 
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which is the desired condition on the The above 

condition is satisfied if rriij = rrii — mj for all i and j. 
[Note, however, that this condition is sufficient but not 
necessary to satisfy Eq. (2.311.] In this case, one can 
redefine the quasienergies via the gauge transformation 


e,- = e,; + TO,; UJ 


(2.32a) 


such that Eq. (2.29) becomes 


R. 

Ri 




(2.32b) 


which yields the desired thermal quasienergy distribu¬ 
tion Pi oc . 


III. ZERO-TEMPERATURE 
NONEQUILIBRIUM STEADY STATES 


In this section, we carry out the program outlined in 
the introduction, and characterize the zero-temperature 
occupations of Floquet bands in a generic class of 
periodically-driven noninteracting systems. We take the 
combined Hamiltonian for the system and the reservoir 
to be given by Eq. (2.1), i.e. 




(3.1a) 


The Hamiltonian H^{t) describes the periodically- 
driven system, which we take for simplicity to be a 
two-band model of noninteracting fermions. (Gen¬ 
eralizing our results to models with more than two 
bands is straightforward.) We will focus on the case 
of monochromatic driving, so that 




(3.1b) 


where A is the driving amplitude and HY){t) depends 
on time via linear combinations of sin(a;t) and cos{u}t). 
The quasienergy bands of the driven system, denoted 
by ej{k) {j = 1,2), are assumed to be gapped (i.e. non¬ 
degenerate for all k) and to host nontrivial topological 
invariants. The system-bath coupling is again taken to 
be in the factorized form of Eq. (2.2), namely 


Hsb=iSB, 


(3.1c) 


We also take the operator S to conserve both momen¬ 
tum and particle number, so that each k mode is ef¬ 
fectively coupled to its own bath, and so hereafter we 
suppress the momentum index k. Finally, Hb describes 
a bosonic bath, with energy eigenstates {v) and eigen¬ 
values El,, which we take to be in equilibrium at zero 
temperature for all time. 

If the Floquet spectrum is nondegenerate, as we as¬ 
sume, then Pij(t) is diagonal at long times, and the 


steady-state occupation probabilities pj := P aji — » oo), 
where pi + P 2 = 1, satisfy the rate equation (2.10), i.e. 


0 — P2 — Ri->-2 Pi- 


(3.2) 


We recall that the transition rates are defined as 

OO 

:= 27r7^ ^ |5'ij(TO)|^g(ej - Ci - tow), 

m— — oo 


Si,{m) := ^ r dte~* (u,(t)| S \uj{t)) . 
T Jo 


(3.3a) 

(3.3b) 


At zero temperature, the w eighted bath density of states 
g{E) defined in Eq. ( |2.5e| ) is given by 

g{E)='£\{0\B\i2)\H{E + E,.). (3.3c) 




We have set the bath ground-state energy Eq = 0, so 
that Ejy^o is strictly positive. For this reason, g{E > 
0 ) vanishes identically at zero temperature, a fact that 
will be of crucial importance below. In our analysis, it 
will be instructive to model g{E) as a power law at low 
energies compared to a very large cutoff scale m- 

g{E)=go\E\^e{-E). (3.4) 


The real exponent 0 < rj < oo classifies the type of bath; 
if ry = I, the bath is referred to as ohmic, while rj > 1 
and ?7 < 1 are referred to as super-ohmic and sub-ohmic, 
respectively. 

It is interesting to note, as po inted out in Sec. |H A[ 
that the rates entering equation (3.2) are invariant un¬ 


der gauge transformations of the form (1.3). While 


these gauge transformations change the ordering of the 
quasienergies, they nevertheless do not change the oc¬ 
cupations of the Floquet states themselves. However, as 
we will see later on, an appropriately engineered reser¬ 
voir is capable of determining an unambiguous ordering 
of the quasienergies and a choice of an effective “Floquet 
ground state.” 

We now turn to an analysis of the transition rates 
(3.3a) that will allow us to determine the extent to which 


a single Floquet state can be populated at zero temper¬ 
ature, given that the system is coupled to a bath. For 
this analysis, it will be convenient to choose a gauge in 
which the ordering of the quasienergies is determined by 
the ordering of the energies of the undriven Hamiltonian 
iJo, in such a way that the separation A := e 2 — ci is 
positive. This gauge can be understood by building up 
the Floquet states perturbatively in A from the eigen- 
basis of Hq m- To do this, we make use of the Fourier 
decomposition of the time-periodic Floquet states, 


E 


girn<^t ^ 


(3.5) 
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and observe that, to zeroth order in this gauge, |m°) 
is nothing but an eigenstate of Hq. (Notice also that 
\uj{nT)) = eigenstates of Hp underly¬ 

ing the Floquet topological insulator.) The nonzero-m 
Fourier components of \uj{t)) arise due to hybridization 
of the spectrum via Hy), and therefore scale with the 
driving amplitude A as 




{XM 


\m\ 


(3.6) 


which follows within perturbation theory to Imj-th order 
in A [32]. Henceforth, we will take A/w to be small either 
on account of a small A or a large lo. In the Appendix, we 
examine the case where A/w is not small, which is much 
less favorable for Floquet topological states. Factoring 
out the A/w scaling from the states |m™), we find that 
one can rewrite the matrix elements in Eq. (3.3b) in the 
following form: 


S,Am)= y (A/w)l'=l+l'=+'"l 

fciy (3.7) 

= (AMI™I s™(A/a;), 

where s'^(X/uj) are regular functions containing the 
terms in the series that depend weakly on A/w or lead 
to a decay of Sij{m) faster than (A/w)!'"! as X/uj —)■ 0. 
With all this in mind, we now analyze the quantity 

P2 _ Ri^2 _ Em l'5'l2(-m)|2g(A-tTOw) 

Pi R2 ^i Eml5i2(+m)Pg(-A + ma;)’ 


which is proportional to the density of excitations above 
the “lowest” quasienergy state in this gauge (ei). If only 
m = 0 above contributes, the excitation density vanishes 
and the lower band is completely filled at long times, as 
it would be at equilibrium, owing to the fact that the 
argument of g{E) in the numerator is positive. While 
model systems that reach such an effective equilibrium 
steady state have been studied |2ll[29H3ll|33|, it is well- 
known that these steady states do not occur for generic 
choices of Hsb and g{E). We ask, instead, whether 
there are any more general mechanisms or limits that 
suppress P 2 lPi- The analysis is simplified if we assume 
that A/w is sufficiently small that we can keep only the 
lowest nontrivial value of \m\ in the sums above. We 
will focus on the case oj > A, since the opposite case 
would likely not yield a topological band structure. 

For uj > A, then only the terms with m < 0 (m < 0) 


contribute to the numerator (denominator) of Eq. (3.8). 


In this case, we find that suppression oip 2 / pi is possib 
in the limit a; ^ A, which yields [c.f. Eq. (3.7)] 


P2 

Pi 


^{X/coW 


’12 


(AMP 


\sl,{X/u;W 




A 


(3.9) 


In addition to the exponent 77 , the behavior of P 2 IP 1 
as a; —)■ 00 depends on the scaling of the quasienergy 
separation A with A and w. We assume the scaling 
A ~ A(A/w)“, for a > 0; for example, the case a = 
1 corresponds to the size of the direct gap predicted 
in graphene coupled to a circularly polarized electric 
field [U |3| . We additionally allow the driving amplitude 
to scale with the frequency, A ~ as it may in some 
physical driven systems Using these scaling 

forms, we find that as a; —)■ 00 the excitation density 

^ - ^-(l-/3)[2-r,(a-Hl)] ^ 

Pi 

SO long as the product 

{1- P)[2-r]{a + l)]>0. (3.10b) 


Noting that the simultaneous requirement of small X/ui 
and large oj restricts /3 < 1, this criterion reduces to 


g < 


2 

a -I- 1 


(3.10c) 


If this condition is not satisfied, then p 2 jpx approaches 
a nonuniversal constant in the high-frequency limit, and 
no single Floquet state is fully populated in the steady 
state. It is important to note that a, and therefore the 
excitation density, is fc-dependent, since the scaling of 
A with A/w varies in momentum space. In particular, 
for k far away from the value at which the minimal 
quasienergy separation occurs, A becomes independent 
of A and w. The w-scaling for this case is obtained by 
setting a = 0 in Eqs. (3.10). 

For a given a and (3 (fixed by the physical realization 
of the system), the above result suggests that the low- 
energy behavior of the function g{E) essentially deter¬ 
mines whether or not a single Floquet band is occupied 
in the limit w 3> A. For example, in the case of graphene 
in a circularly-polarized electric field (a = 1 , /3 = 0 ), 
there is a critical value r^c = 1 (i.e., ohmic dissipation) 
that separates the power-law decay of P 2 IP 1 from the 
aforementioned nonuniversal behavior. Therefore, an 
ohmic bath already violates Eq. (3.10c) for graphene in 
a circularly-polarized electric field, and the population 
of the bands near the AT-point is not controllable by 
increasing ui if the bath is ohmic. 

In cases where the excitation density decays as a 
power law at large freqnencies, one must still take care 
to determine whether the resnlting steady state has the 
desired characteristics of the topological Floquet bands. 
As UJ increases, the Floquet effective Hamiltonian Ely 
may approach Hq as 1/uj or faster. If the power-law 
decay of P 2 IP 1 is faster than this approach, then the 
suppression of excitations can still occur in a regime 
where the Floquet bands are topological. The situation 
can be improved by allowing a scaling of the amplitude 
A ~ for /3 7 ^ 0 [13], but the value of the exponent 
















E 



FIG. 3: (Color online) Suppressing unwanted transitions be¬ 
tween Floquet states by design of the function g{E). De¬ 
pleting the density of states in the neighborhood of such a 
transition ensures that there is no corresponding transition 
in the bath to compensate, essentially forbidding it. 


j3 must be balanced against an appropriate value of rj 


[c.f. Eq. (3.10c I] in order for the desired suppression 
to take place. Furthermore, it is important to keep in 
mind that, depending on the exponents a, j5, and 77 , the 
power law decay of P 2 IP 1 can be very slow, so that a 
finite density of excitations remains even at high fre¬ 
quencies compared to all system energy scales. In order 
for excitations to be completely suppressed, one must 
have 


Teff < A, (3.11a) 


where the frequency- and momentum-dependent effec¬ 
tive temperature Tea is defined as 




A 

In(pi/P2)' 


(3.11b) 


Note that this effective temperature arises even though 
the bath itself is at zero temperature, and can therefore 
be understood as a signature of heating effects due to 
the driving. 


IV. THE ROLE OF BATH ENGINEERING 


The sensitive dependence on the exponent rj of the 
large-frequency scaling of the excitation density already 
demonstrates the crucial role that the bath plays in sta¬ 
bilizing Floquet topological states in open systems. We 
close by commenting on two additional ways in which 
the bath and its coupling to the system can be engi¬ 
neered in order to further favor the suppression of P 2 /P 1 • 
One way to suppress excitations is to engineer the 
spectrum of the bath i tself so that the function g{E) 
appearing in Eq. ( 3.3a | does not have a simple power- 
law form as in Eq. (3.4), but instead drops to zero in 


a neighborhood oi E = Aq — w, where Aq := A(fco) 
is some reference value of the quasienergy separation. 


as in Fig. (Such a scenario could be envisioned if, 
for example, the bath consists of quantized electromag¬ 
netic radiation in a cavity, whose size could be tuned 
to achieve the desired effect.) If the width of the dip 
in g{E) is on the order of the width of the upper band, 
then excitations can be suppressed even if w is not much 
larger than Aq. Indeed, if w > Aq and X/ui is sufficiently 
small that only the first few terms in the sums over m 
in Eq. ( |3.8[ ) are kept, then the suppression of g{E) near 
this value eliminates the terms with m 7^ 0 up to or¬ 
der (A/w)^. However, it is important to point out that 
the excitation density in this case still exhibits at best 
power-law decay at la rge fr equencies, with appropriate 
modifications to Eqs. (3.10) arising from keeping terms 


other than \m\ = 0,1 in Eq. ( |3.8[ ). To completely eradi¬ 
cate excitations to all orders in A/w, one must engineer 
dips at energies E = Aq — muj for all m > 0. 

Of course, by placing the dip at E = — Aq — w, one 
can also use this mechanism to populate what we have 
referred to as the upper band in this choice of gauge. 
While this scenario looks like a population inversion, 
one can of course perform an appropriate transforma¬ 
tion of the form (1.3) to reorder the Floquet bands in 
such a way that the populated band is the lowest. This 
example indicates that, in certain cases, the reservoir 
can “choose” a preferred gauge in which the system ap¬ 
pears to be (nearly) at equilibrium. 

The system-bath coupling is another quantity that 
could potentially be manipulated in order to suppress 
excitations. Indeed, certain system-bath couplings are 
known to yield relaxation to steady states that feature 
filled Floquet bands. For example, if Hsb is chosen 
in s uch a way that P''' {t)HsBP{t), where P{t) satisfies 
Eq. (1.2) with H{t) = Hs{t), is time-independent, then 
the system described by the total Hamiltonian E[ (t) de¬ 
fined in Eq. (|2.1|) reaches an effective thermal equilib¬ 


rium with respect to the eigenvalues of Pip (see Sec. H B 
and Refs. 


Given sufficient control over the 
system-bath coupling, one could attempt to engineer 
such a situation, at least to some order in A/w, by de¬ 
signing an iJsB that, say, commutes with the lowest non¬ 
trivial Fourier harmonic of P(t), or even by engineering 
an appropriate time dependence in i?sB to cancel the 
time-dependence in P{t) to some order. 


V. SUMMARY AND CONCLUSION 


We have argued in this work that the possibility of 
stabilizing a Floquet topological state with a low den¬ 
sity of excitations is heavily constrained by the coupling 
to a thermal reservoir. Using scaling arguments, we 
demonstrated that, even in the limit of weak driving 
and/or high frequenscy, the bath density of states has 
tremendous influence on whether or not excitations are 
suppressed as w —)■ 00 . We also suggested ways of de- 
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signing the bath and its coupling to the system in order 
to suppress excitations. 

Our results suggest that it is at best difficult, and 
at worst impossible, to engineer a periodically-driven 
quantum system whose steady state resembles the zero- 
temperature ground state of some target topological 
phase. However, even out of equilibrium, there is reason 
to believe that nontrivial features, such as topological 
indices |3S], edge states [S], and (approximately) quan¬ 
tized transport HUH], survive in both isolated and open 
systems. Indeed, there is already experimental evidence 
to this effect in cold atomic gases [331 I3S]- We em¬ 
phasize, however, that it is precisely in the deviations 
from the resemblance to equilibrium systems where the 
newest physics lies. For example, interacting versions 
of these models [B] could be used as platforms to probe 
fractionalized excitations out of equilibrium. 
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negligible. Indeed, as A/w —)■ oo, the Floquet states can 
become chaotic, so that the Sij{m) may be regarded 
as essentially random variables, whose magnitudes need 
not decay quickly as \m\ becomes large. For this reason, 
the sums in the numerator and denominator of Eq. (3.81 
generically diverge in the limit X/ui —>■ oo, and the ratio 
of transition rates is indeterminate. One can, however, 
identify constraints on the amplitudes |S'i 2 (to)P such 
that the ratio converges to a definite finite value. In 
particular, if 


|5'i2(±|m|)p < (const.) x ^^+^,+5 (^1) 


as \m\ —>■ oo for any positive real number <5, then both 
series are bounded from above by a convergent series, 
and therefore the ratio has a definite value. This is true 
even for infinitesimally small (5 —> 0+. 

Even if the sums in the numerator and denomina¬ 


tor are divergent, the ratio (3.8) can approach a finite 


value for system-bath coupling operators S such that 
\Sij{m)\^ = |5'ij(—m)p, due to a symmetry. To see 
this, let us drop the m = 0 term in the denominator 
and rewrite Eq. (|3.8|) for w > A as 


^ _ Em <0 \Si 2 {\m\)\^ {\m\uj - A)^ 

Pi Em<ol‘5'l2(|w|)p(k|w + A)''" 


(A2) 


Note added — During preparation of this manuscript, 
we became aware of Ref. [25l which also discusses the 
possibility of stabilizing Floquet topological states with 
couplings to particular appropriately-engineered baths. 

Appendix: Limit of strong driving 

In the case of strong driving (A/w » I), there 
are many values of m for which Sij{m) can be non- 


For large |to|, the summands in the numerator and de¬ 
nominator become identical. Therefore, if |S'i2(|»7i|)P 
is finite for sufficiently large |m|, the ratio of the two 
sums approaches I from below as \m\ grows. When this 
occurs, the system approaches an infinite effective tem¬ 
perature — all Floquet states are occupied with equal 
probabilities, despite the fact that the bath is held at 
zero temperature. If there exists some Irumaxl such that 
|5'i2(|m.max|)P = 0, then the ratio takes on a finite value 
that is bounded from above by unity. 
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